Long-range ferromagnetic dipolar ordering of high-spin molecular clusters 
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(Dated: February 6, 2008) 

We report the first example of a transition to long-range magnetic order in a purely dipolarly inter- 
acting molecular magnet. For the magnetic cluster compound Mn604Br4(Et2dbm)e, the anisotropy 
experienced by the total spin S = 12 of each cluster is so small that spin-lattice relaxation remains 
fast down to the lowest temperatures, thus enabling dipolar order to occur within experimental 
times at T c = 0.16 K. In high magnetic fields, the relaxation rate becomes drastically reduced and 
the interplay between nuclear- and electron-spin lattice relaxation is revealed. 
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Few examples of long-range magnetic order induced 
by purely dipolar interactions are known as yet |Jx|, ||. 
Therefore, the possibility to study such phase transitions 
and the associated long-time relaxation phenomena in de- 
tail in high-spin molecular cluster compounds, with vary- 
ing crystalline packing symmetries and different types of 
anisotropy, presents an attractive subject j|. However, 
for the most extensively studied molecular clusters so- 
far such as Mni2, Fe 8 and Mn 4 j|, |[ ||, J?J the uniaxial 
anisotropy experienced by the cluster spins is very strong. 
Consequently, the electronic spin-lattice relaxation time 
Tf becomes very long at low temperatures and the clus- 
ter spins become frozen at temperatures of the order of 
1 K, i.e. much higher than the ordering temperatures 
T c ~ 0.1 K expected on basis of the intercluster dipo- 
lar couplings ||. Although quantum tunneling of these 
cluster spins has been observed [|, ^j, [?J , and could in 
principle provide a relaxation path towards the magnet- 
ically ordered equilibrium state ||, the associated rates 
in zero field are extremely small (< 100 Hz). For these 
systems, tunneling only becomes effective when strong 
transverse fields B t are applied to increase the tunneling 
rate. Although the tunability of this rate and thus of Tf 2 
by B t could recently be demonstrated for Mni 2 , Fes and 
Mn 4 0, no ordering has yet been observed. 

The obvious way to obtain a dipolar molecular mag- 
net is thus to look for a high-spin molecule having suf- 
ficiently weak magnetic anisotropy and negligible inter- 
cluster super-exchange interactions. Here we report data 
for Mn604Br4(Et2dbm)6, hereafter abbreviated as Mng 
P). The molecular core of Mng is a highly symmetric 



octahedron of Mn 3+ ions (with spin s — 2) that are 
ferromagnetically coupled via strong miro-cluster super- 
exchange interactions (Fig. 1). Accordingly, the ground 
state is a S = 12 multiplet and the energy of the near- 
est excited state is approximately 150 K higher ||. The 
unit cell is monoclinic, with space group Pc and con- 
tains 4 molecules || bound together only by Van dcr 
Waals forces, /nier-cluster superexchange is therefore 
negligible and only dipolar interactions couple the clus- 
ter spins. The net magnetocrystalline anisotropy of this 
cluster proves to be sufficiently small to enable measure- 
ments of its equilibrium magnetic susceptibility and spe- 
cific heat down to our lowest temperatures (15 mK). 




FIG. 1: (a) ORTEP representation {J of the structure of 
Mno04Br4(Et2dbm)6; (b) detail of the symmetric octahe- 
dral core, containing six ferromagnetically coupled Mn 3+ ions, 
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The main results of this paper follow. We have ob- 
served that the magnetic clusters do undergo a transi- 
tion to a long-range ferromagnetically ordered state at 
T c = 0.161(2) K. This transition can be observed within 
experimental times (1 — 100 s) because the magnetic re- 
laxation is not too slow at very low temperatures. We 
have obtained further information on the magnetic re- 
laxation by applying a magnetic field. It turns out that 
the relaxation rate vanishes exponentially as the field in- 
creases. The interplay between nuclear- and electron-spin 
lattice relaxation is also revealed. 

Polycrystalline samples of Mng were prepared as in 
Ref. . The specific heat of a few milligrams of sam- 
ple, mixed with Apiezon grease, was measured at low-T 
in a home-made calorimeter M that makes use of the 
thermal relaxation method. An important advantage of 
this method is that the characteristic time r e of the ex- 
periment (typically r e ~ 1 — 100 seconds at low-T) can 
be varied by changing the dimensions (and therefore the 
thermal resistance) of the Au wire that acts as a thermal 
link between the calorimeter and the mixing chamber of 
the 3 He- 4 He dilution refrigerator. The ac susceptibility 
data were taken between 0.015 and 4 K with a mutual 
inductance bridge. The frequency was varied from 230 to 
7700 Hz. Magnetic data for T > 1.8 K were taken with 
a SQUID magnetometer. 

We first discuss the magnitude of the anisotropy. The 
spin Hamiltonian for the molecule is: 

H= -DS 2 Z - gp B B ■ S. (1) 

Previously published magnetic data for T > 1.8 K al- 
ready indicated an upper limit of about 0.01 K for D/ks 
In order to obtain an independent estimate of D, 
high-frequency ESR data were taken in the range 95—380 
GHz. Owing to the combination of very small D and 
large S, as well as the presence of a signal at g = 2.00 
arising from a minute amount of Mn 2+ impurity (often 
seen in ESR of Mn 3+ compounds), the interpretation of 
the spectra was not fully conclusive. Nevertheless, sig- 
nals with a clearly visible structure on the low-field end 
of the spectra could be obtained. It could be identified as 
fine structure originating from zero- field splitting (ZFS), 
since it was independent of field and frequency Simula- 
tions of the spectra performed using Eq. (|l|) agree well 
with the experiment taking \D\/ks ~ 0.03 or ~ 0.05 K, 
depending on the sign of D (which could not be unequiv- 
ocally determined). Although a smaller rhombic compo- 
nent could be present, the data do not justify a more 
elaborate fitting. 

The isotropic character of the molecular spin might 
seem paradoxical at first, considering that the individ- 
ual Mn 3+ ions, being Jahn- Teller ions, experience strong 
anisotropy. However, the spin Hamiltonian for the clus- 
ter is determined by the vectorial addition of the local 
tensors of the individual atoms, which can give rise to a 
low net anisotropy for highly symmetric molecules such 
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FIG. 2: Ac susceptibility of Mri6 at zero field and four differ- 
ent frequencies. Data measured above 1.8 K in a commercial 
SQUID magnetometer are also shown (x). Inset: x' corrected 
for demagnetization effects. The full line shows the best fit 
of a Curie- Weiss law to the data for T > 0.3 K, yielding 
= 0.20(3) K. 

as Mng (cf. Fig. 1), no matter how large the ZFS of 
the constituting atoms ]1C[| . In fact, the possibility to 
tune the net anisotropy of the cluster spin is one of the 
attractive properties of molecular superparamagnets. 

As first evidence for ordering of the magnetic moments 
we show zero-field ac susceptibility data in Fig. 2. The 
real part \' shows a sharp maximum at T c — 0.161(2) 
K. We found that x' at T c is close to the estimated limit 
for a ferromagnetic powder sample, l/(pN s + p sa mN c ) ~ 
0.14±0.02 emu/g, where p = 1.45 g/cm 3 and p sam ~ 0.45 
g/cm 3 are respectively the densities of bulk Mng and of 
the powder sample, N s — Att/3 is the demagnetizing 
factor of a crystallite, approximated by a sphere, and 
N c ~ 2.51 is the demagnetizing factor of the sample 
holder. This indicates that Mng is ferromagnetically or- 
dered below T c . Susceptibility data \i corrected for the 
demagnetizing field (xi = x'/I 1 - (pN 8 + p sam N c )x']) fol- 
low the Curie- Weiss law x = C/(T — 8) down to approxi- 
mately 0.3 K, with C = 0.034(1) cm 3 K/g and 8 = 0.20(3) 
K. The constant C equals, within the experimental er- 
rors, the theoretical value for randomly oriented crystals 
with Ising-like anisotropy Na9 2 P%S(S + l)/3fcsP m = 
0.0332 cm 3 K/g, where S = 12, g = 2, and the molec- 
ular weight P m — 2347.06. The positive 8 confirms 
the ferromagnetic nature of the ordered phase. From 
the mean-field equation 8 = 2zJ e ffS(S + l)/3fcs, we 
estimate the effective inter-cluster magnetic interaction 
J e ff w 1.6 x 10~ 4 K, and the associated effective field 
H e ff = 2zJ e f fS / gpB = 3.5 x 10 2 Oe coming from the 
z = 12 nearest neighbors. 

The maximum value of x' is seen to vary only weakly 
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FIG. 3: Electronic specific heat data (O) of Mri6 at zero 
field. The dotted line is the Schottky anomaly calculated 
with D = 0.013 K. The full line is the Monte Carlo (MC) 
calculation for an orthorhombic lattice of 1024 Ising spins with 
periodic boundary conditions. For each point, we performed 
2 x 10 4 MC steps per spin. 

with cj, which we attribute to the anisotropy. The total 
activation energy of Mng amounts to DS 2 w 1.5 K, i.e., it 
is about 45 times smaller than for Mni2 . Accordingly, one 
expects the superparamagnetic blocking of the Mn@ spins 
to occur when T ~ TB(Mni 2 )/45, that is below ~ 0.12 K. 
In other words, for T — > T c , the approach to equilibrium 
begins to be hindered by the anisotropy of the individual 
molecular spins. We stress however, that the frequency 
dependence of x' observed here is very different from that 
of the well known anisotropic superparamagnetic clusters 
or that of spin glasses. Below T c , \' decreases rapidly, 
as expected for an anisotropic ferromagnet in which the 
domain-wall motions become progressively pinned. The 
associated domain-wall losses should then lead to a fre- 
quency dependent maximum around T c in the imaginary 
part, x"t as seen experimentally. Indeed, although the 
Mn6 spins can be considered as nearly isotropic at high 
temperatures, the anisotropy energy (~ 2DS) is of the 
same order as the dipolar interaction energy /i 2 /r 3 ~ 0.1 
K between nearest neighbor molecules. Thus the order- 
ing should be that of an Ising dipolar ferromagnet. 

Additional evidence for the ferromagnetic transition 
at T c is provided by the electronic specific heat c e [ pd| , 
shown in Fig. 3, which reveals a sharp peak at 0.15(2) K. 
Numerical integration of c e /T between 0.08 K and 4 K 
gives a total entropy change of about 3Ak B per molecule, 
very close to the value for a fully-split 5=12 spin multi- 
plet (i.e., k B ln(25"+l) = 3.22/c B ). We thus may attribute 
the peak to the long-range order of the molecular spins. 
We note, however, that at T c the entropy amounts to 
about lk B per spin, showing that only the lowest energy 
spin states take part in the magnetic ordering. 

We have performed Monte-Carlo (MC) simulations for 
an S — 12 Ising model of magnetic dipoles on an or- 



thorhombic lattice with axes a x = 15.7 A, a y = 23.33 A, 
and a z = 16.7 A, which approximates the crystal struc- 
ture of Mn6 . The model includes dipolar interactions as 
well as the anisotropy term —DS 2 given in Eq. ([!]). Ne- 
glecting interactions, the ZFS of the S = 12 multiplct 
produced by this crystal field term lead to a Schottky 
anomaly in c e , as shown by the dotted curve in Fig. 3. 
The fit in the range above 0.5 K yields D/k B — 0.013 
K, in good agreement with our above estimates. The 
intermolecular dipolar interactions remove the remain- 
ing degeneracy of the | ± m) spin doublets. The MC 
simulations show that the ground state is ferromagnet- 
ically ordered, as observed, and predict a shape for c e 
that is in very good agreement with the experiment. In 
Fig. 3, we show c e calculated assuming all molecular 
easy (z) axes to point along a z , i.e. one of the two nearly 
equivalent short axes of the actual lattice. Similar re- 
sults were obtained for other orientations chosen for the 
anisotropy (z) axis. We note that the Ising simulations 
give T c — 0.22 K, which is slightly higher than the exper- 
imental T c — 0.161(2) K. This difference may be related 
to the finite value of the anisotropy. Model calculations 
for this crystal structure, assuming classical anisotropic 
Heisenberg spins with varying anisotropy jl2| , show that 
different ferromagnetic ground states are possible, de- 
pending on the competition between local crystal field 
effects and long-range dipolar interactions. The varia- 
tion of T c with anisotropy, as well as the form of the 
calculated and observed specific heat anomaly, are spe- 
cific for dipolar interactions, and differ widely from the 
analogues for the usual superexchange ferromagnets |l3[ . 

We next turn to the specific heat data obtained in 
varying magnetic field B, plotted in Fig. 4. Even for 
the lowest B value, the ordering anomaly is fully sup- 
pressed, as expected for a ferromagnet Jl^]. Accordingly, 
we may account for these data with the Hamiltonian (1) 
neglecting dipolar interactions. The Zeeman term splits 
the otherwise degenerate | ±m) doublets, and already for 
B ~ 0.5 T the level splittings become predominantly de- 
termined by B, so that the anisotropy term can also be 
neglected. As seen in Fig. 4, the calculations performed 
with D = reproduce the data quite satisfactorily at 
higher temperatures (dotted curves). 

However, when the maxima of the Schottky anomalies 
are shifted to higher T by increasing B, an additional 
contribution at low T is revealed. It is most clearly vis- 
ible in the curves for 1 T< B < 2.5 T, and varies with 
temperature as cT 2 /R = 4 x 10~ 3 . We attribute this to 
the high-temperature tail of the specific heat contribu- 
tion Cnud arising from the Mn nuclear spins (/ = 5/2), 
whose energy levels are split by the hyperfine interaction 
with the Mn 3+ electronic spins s. This interaction can be 
approximated by Hhf — sAhftni, where Ahf is the hy- 
perfine constant and m/ is the projection of the nuclear 
spin along the electronic spin. At high temperatures (i.e., 
when A hf s < k B T) c nucl /R ~ \A 2 hf s 2 I{I + l)T- 2 @. 



4 



t 1 1 — i — i i i i I 1 1 1 — i — i i i i I 1 1 1 — r 




T(K) 

FIG. 4: Magnetic specific heat of Mne for various applied 
fields. Dotted and dashed lines: calculated equilibrium spe- 
cific heats of the electronic and the nuclear spins, respec- 
tively. Full lines: time-dependent total specific heat calcu- 
lated taking into account the nuclear spin-lattice relaxation. 
Inset: magnetic field Bb needed to take the nuclear spins off- 
equilibrium, o, obtained from c nuc i vs B isotherms; •, from 
Cnuci vs T at constant field. 

Taking A^f = 7.6 mK as used previously to simulate 
ESR spectra measured on a Mn 4 cluster pl|, we obtain 
the dashed line of Fig. 4. This contribution was sub- 
tracted from the zero-field data shown in Fig. 3. 

A remarkable feature of the experimental data that is 
not reproduced by these calculations is that, at the lowest 
T, the nuclear specific heat drops abruptly to about 10 -5 
J/gK. The temperature Tb where the drop occurs de- 
pends on B but also on the characteristic time constant r e 
of our (time-dependent) specific heat experiment: the de- 
viation from the (calculated) equilibrium specific heat is 
found at a lower T when the system is given more time to 
relax. We conclude that the drop indicates that nuclear 
spins can no longer reach thermal equilibrium within time 
r e . We may write c„ uc ;(r e ) = c** rf [l - exp(-T e /Ti)], 
showing that the transition should occur when the nu- 
clear spin-lattice relaxation time T\ becomes of the or- 
der of T e . These transitions to non-equilibrium provide 
therefore direct information on the temperature and field 
dependence of T\, which can be related to the fluctu- 
ation of the transverse hyperfine field, as produced by 
the phonon-induced transitions between different levels 
of the electronic spin |ll| . 

At low T and high B, only the ground and the first ex- 
cited states, m — +12 and m = +11, need be considered. 
We may write l/Ti ~ (1/to) exp{ — [g^sB]/kBT}, where 
l /to plays the role of an attempt frequency for the elec- 
tron spin transitions. Clearly, the nuclear spins can be 
taken out of of equilibrium either by decreasing T down 
to Tb at constant field (as in Fig. 4) or by increasing 
B up to a given value Bb at constant T. This is indeed 



observed experimentally (not shown). The effect of the 
field is just to polarize the electronic spins, which reduces 
the fluctuations of the hyperfine field, thus effectively dis- 
connecting nuclear spins from the lattice. For a given r e , 
Bb must increase linearly with Tb, which is confirmed by 
the experimental data plotted in the inset of Fig. 4. The 
slope gives r ~ 3 x 10 -4 s. Using this value, we have cal- 
culated the time-dependent c nuc i, shown as the full lines 
in Fig. 4, and seen to be in reasonable agreement with 
the experimental data at all T and B. 
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